Magnetism in Exact Exchange Density Functional Theory 
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The magnetic properties of the intermetallic compound FeAl are investigated using exact exchange 
density functional theory. This is implemented within a state of the art all-electron full potential 
method. We find that FeAl is magnetic with a moment of 0.70 ^b, close to the LSD A result of 
0.69 fj,B. A comparison with the non-magnetic density of states with experimental negative binding 
energy result shows a much better agreement than any previous calculations. We attribute this to 
the fine details of the exchange field, in particular its asymmetry, which is captured very well with 
the orbital dependent exchange potential. 

PACS numbers: 75.10.-b, 71.10.-w 



INTRODUCTION 

The intermetallic compound FeAl presents an inter- 
esting system where one has both a complex interplay 
between structure and magnetism, and a possible fail- 
ure of the Local Spin Density Approximation (LSDA) to 
adequately account for the electronic structure [l], y, |3| . 
Disentangling these two strands poses a severe challenge 
for theory. The controversy surrounding this material 
turns on the issue of whether completely ordered FeAl is 
a magnetic or non-magnetic intermetallic, and whether 
this question is at all relevant to the experimental pic- 
ture. 

Experimentally, FeAl obeys an extended Curie- Weiss 
law {x~Xo)~^ = C{T—<d) down to low temperatures, the 
signature of a paramagnet with finite (local) moments. 
Mossbauer studies Q show that it is only the antisite Fe 
atoms (Fe atoms on the Al sublattice) which contribute 
to this behaviour, the Fe sublattice being non-magnetic. 
The rather weak ordering in FeAl means both that a 
range of defects are thermodynamically stable, and that 
any crystal of FeAl will have some concentration of these 
defects. Thus it appears that the question of whether 
completely ordered FeAl is magnetic or not cannot be 
answered by experiment. 

Nevertheless, the question of the magnetic state of 
completefy ordered FeAl has attracted a lot of atten- 
tion 0, U 0) S IE ■ This is based on the speculation 
that it should be non-magnetic and the fact that 



LSDA calculations find it to be magnetic with a moment 
of around 0.69 fis- The foundation of this speculation 
appears to be simply that since the magnetism is carried 
by defects, taking these defects out will leave behind the 
non-magnetic FeAl sublattice, hence non-magnetic FeAl. 
However, this excludes that an Fe antisite-Fe sublattice 
interaction may be responsible for making the Fe sublat- 
tice non-magnetic. 

The first attempt to go beyond the LSDA was un- 
dertaken by Mohn et. al who employed the LDA-I-U 
method. The usefulness of this approach was suggested 
by the fact that both the Cg and t2g bands are rather 
flat, and hence an improved scheme for the correlation 
is needed. Although these authors found a non-magnetic 
solution, a more recent LDA-I-U calculation 1] disagrees 
with this result and finds, again, FeAl to be magnetic. 
The latter work instead suggested that the magnetism 
may be suppressed by critical spin fluctuations, as is 
seen for example in the intermetallic compound NisGa 
7;]. The dynamical version of LDA-I-U, the dynamical 
mean field theory (DMFT)0], does include some spin- 
fluctuations and using this it was found that FcAl is non- 
magnetic. 

In complete contrast to these approaches Smirnov et. 
al instead attempted a more realistic modeling of 
the experimental defect structure within the LSDA. To 
achieve this they employed the Coherent Potential Ap- 
proximation (CPA), which allows for the treatment of 
both structural and magnetic disorder. With an Fc anti- 
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site concentration similar to that deduced experimentally 
they found that a spin disordered local moment (DLM) 
solution is then energetically favoured. Crucially, in this 
state the Fe sublattice is non-magnetic, with the mag- 
netism carried only by the antisite atoms, exactly as seen 
in experiment. 

Although it is clear that the complex defect structure 
must be taken into account in order to model correctly 
the experimental reality, one can note that it is not clear 
that the magnetism of FeAl is correctly described by the 
DLM state. Evidence for this comes from the fact that 
Fe clusters in FeAl order with a non-collinear structure. 
Thus the Fe antisite-Fe sublattice interaction is not fer- 
romagnetic, and hence the DLM state should be lower in 
energy. This does not, however, rule out the possibility 
that some complex non-collinear structure is even lower 
in energy, and that in this state the Fe lattice is again 
magnetic. 

In this situation, one might wish for a more experi- 
mentally rigorous comparison than simply whether FeAl 
is non-magnetic or not. Unfortunately, and perhaps sur- 
prisingly, there is a lack of experimental spectroscopic 
data for this system. To the best of our knowledge, only 
one experimental article on this aspect of FeAl is to be 
found in the literature 0. 

In this work we deploy the EXact eXchange 
(EXX) method, implemented within an all-electron full- 
potential code EXITING jl], to study both the mag- 
netism and spectroscopic properties of FeAl. 

METHOD OF CALCULATION 

Recently we have developed EXX equations suitable 
for non-collinear magnets [9|. Within this formalism the 
exchange term is treated exactly within Kohn-Sham den- 
sity functional theory (KS-DFT). The resulting exchange 
potential is self-interaction free and depends explicitly on 
the KS orbitals, which are two component Pauli spinors. 
This formalism leads to coupled integral equations for 
the exchange potential and magnetic field. In order to 
avoid an expensive set of fixed spin moment calculations 
this method is needed even for the study of coUinear 
systems 0- In the present work we use an all elec- 
tron full-potential linearized augmented plane- wave (FP- 
LAPW) method for the self-consistent EXX calculations. 
In this scheme, the unit cell is divided into non overlap- 
ping muffin-tin spheres around the atomic nuclei and a 
remaining interstitial region. The corresponding wave 
functions are expressed as linear combinations of atomic 
like functions inside the spheres and plane waves in the 
interstitial region. The calculations used a RKmax = 7 
cutoff criterion for the number of plane waves in the ba- 
sis. The wave function is expanded up to ^ = 10 and 
the potentials up to Z = 6. All calculations are done for 
the experimental lattice parameter of 5.496 au. We have 



TABLE I: Magnetic moment of FeAl in Bohr magneton. The 
FP-EXX and FP-LSDA are our results. The LDA+U results 
from Refs and Q are for U=5eV. The DMFT results are 
from Ref. 1 and use U=2eV. 

Expt. FP-EXX FP-LSDA LDA+Uf4] LDA-fUfl] DMFT 
? 070 069 070 



used a mesh of 64 k-points in the irreducible part of the 
Brillouin zone. 



RESULTS 

The EXX calculations find a magnetic ground state 
with a moment of 0.70 /z_b, very close to the LSDA value 
0.69/i_B. It would be interesting at this point to compare 
the density of states (DOS) for FeAl with the experi- 
ments. A comparison between the magnetic DOS and 
the experimental negative binding energy results of Sha- 
banova et. al will not be meaningful, since, as has been 
discussed above, the Fe sublattice is non-magnetic in ex- 
periment. However, since only 2%-3% of the Fe atoms 
in the FeAl sample will actually be high moment antisite 
defect atoms, their contribution to the binding energy 
signal will be quite small. Hence a rather good approxi- 
mation to model the experimental results will be to sim- 
ply make the comparison with non-magnetic FeAl. In 
principle, of course, one would wish to include the com- 
plex defect structure within the EXX approach, however, 
since this requires either very large supercells or the use 
of the CPA, this is not feasible at present. 

In the lower panel of Fig. 1 we show the results of an 
EXX calculation of the DOS. A most remarkable feature 
can immediately be seen which is the shift downwards of 
the non-bonding t2g states by almost 1.5 eV. This is quite 
different from the results obtained using the LDA-I-U 
and the DMFT 0] methods. Like the LDA+U method 
we see a broadening of the d-band, in contrast to the 
DMFT density of states, where the d-band is narrowed 
by the mass renormalization effects. This spread of re- 
sults by different methods makes a comparison with ex- 
periments most interesting. In the upper panel of Fig. 
1 we show the LDA and EXX DOS, as well as that re- 
ported in Refs.iandll Note that for ease of comparison 
we have convolved our DOS with a Gaussian of width 0.5 
eV. Clearly, the peak of the experimental DOS and that 
given by EXX are in much better agreement than any 
other method. A further interesting feature of the EXX 
DOS is a pronounced anti-bonding eg peak at around 4 
eV, a feature that should be detectable in inverse pho- 
toemission experiments. 

This enhanced splitting of the t2g and e^ manifolds over 
LSDA values is the result of the asymmetry in the EXX 
potential. This is demonstrated in Fig. 2 where the EXX 



3 




FIG. 1: Total density of states for FeAl. The LDA results are 
shown in red, EXX results with only I = part of potential in 
blue, and full potential EXX in black. The t2g (gray) and e^ 
(brown) manifolds are shown with solid shading for EXX and 
hatched shading for LDA. At the top the green curve is the 
experimental negative binding energy from Ref. 5] alongside 
the theoretical Gaussian convolved DOS. The cyan curve is 
the data extracted from Ref. and orange from Ref. 



exchange potential (black) along with the LDA exchange 
potential (red) are displayed. A qualitative difference be- 
tween the two can be seen in that the EXX potential is 
strongly asymmetric around the Fe atoms. Looking at 
the potential along an Al-Fe-Fe path, indicated by the 
arrows in Fig. 2, one notes structures reflecting the exis- 
tence of atomic shells. These features are a consequence 
of the natural orbital dependence of the EXX potential, 
but are completely absent in the LDA potential which 
is constructed from only local knowledge of the density. 
Dips in the EXX potential will mean an effective attrac- 
tive force which will act to cause localization, the pro- 
nounced dip in the Fe-Fe direction, for example, is most 
likely responsible for the sharp eg peak in the EXX DOS. 

The importance of this asymmetry in the EXX poten- 
tial may be gauged by retaining only the spherical part 
in the muffin-tin spheres during a self-consistent calcula- 
tion. The resulting EXX potential (blue line. Fig. 2 then 
lacks the asymmetry of the full EXX potential, and the 
resulting DOS is much closer to that of the LDA, and in 
disagreement with experiment. 

This agreement of the EXX DOS with experiment can 
be considered somewhat surprising in that previous EXX 
calculations have not found this to be so for the transi- 
tion metals Fe, Co, and Ni. Here a quite striking dis- 
agreement with experiment, as well as the LSDA, was 
noted Although these calculations used the atomic 
sphere approximation, as well as an inaccurate formula- 
tion of EXX for magnetism, our results are qualitatively 



FIG. 2: Upper image shows contour plots of the EXX (black) 
and LDA (red) exchange potential in the [110] plane over two 
unit cells with contour spacing 3.6 eV. In the lower graph, 
the same potentials are plotted along the line indicated in the 
upper graph, with the blue curve corresponding to the I = 
EXX. An offset is added to the EXX potentials to bring them 
into coincidence with the LDA potential. 
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similar dim. One finds the splitting of the t2g and Cg 
states for the partially occupied spin down channel is far 
to large, with the occupied spin up channel being pulled 
up to 3 eV below the Fermi level. One can speculate 
upon the reason for the good agreement between EXX 
and experimental DOS found in this work. Since one of 
the main effects of the missing correlation will be simply 
to screen the exchange field, one can note that the pres- 
ence of Al atoms in this material leads to an effective 
screening by the delocalized sp Al states. The validity 
of this proposal will be investigated in a future work on 
screened EXX. 

In conclusion, we find that FeAl prefers to be ferromag- 
netic in both EXX and LSDA calculations. However, the 
non-magnetic EXX DOS is in very good agreement with 
experiment, while the non-magnetic LSDA DOS is not. 
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Thus are calculations are in agreement with the sugges- 
tion of Smirnov et al. that the non-magnetic state of the 
Fe sublattice results from the complex defect structure. 
However, it is also clear that the LSDA exchange and 
correlation is not sufficient to reproduce (or understand) 
all the properties of this interesting system. 
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